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ABSTRACT 

Non-thermal TeV 7 -ray emission within a multiparsec has been observed from 
the center region of our Galaxy. We argue that these 7 -rays are the result of tran¬ 
sient activity of the massive black hole Sgr A* that resides at the Galactic Genter. 
Several thousand years ago, the black hole may have experienced an active phase 
by capturing a red giant star and forming an accretion disk, temporarily be¬ 
having like an active galactic nucleus. A powerful jet, which contains plenty of 
high speed protons, was launched during the process. These runaway protons 
interact with the dense ambient medium, producing TeV 7 -ray emission through 
the 7 r°-decay process. We show that the total energy deposited in this way is 
large enough to account for observations. The diffusion length of protons is also 
consistent with the observed size of the TeV source. 

Subject headings: Galaxy: center - Galaxies: active - Galaxies: jets - Galaxies: 
gamma rays - Accretion disk - Black holes 


1. Introduction 

It is known that there are many remarkable high-energy sources harbored in the Galactic 
center (GG) region (Melia & Falcke 2001). Recently, TeV 7 -ray emission from the direction 
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of the GC has been reported by three independent groups, Whipple (Kosack et ah 2004), 
CANGAROO (Gollaboration of Australia and Nippon for a Gamma Ray Observatory in the 
Outback; Tsuchiya et ah (2004)), and HESS (High Energy Stereoscopic System; Aharonian 
et ah (2004)). At least four potential candidates are suggested for this TeV y-ray emission, 
which are the black hole Sgr A* (Aharonian & Neronov 2005a; Atoyan & Dermer 2004; 
Levinson 2000), the compact and powerful young supernova remnant (SNR) Sgr A East 
(Grocker et ah 2005), the dark matter halo (Horns 2005; Profumo 2005; Gnedin & Primack 
2004; Ellis et ah 2002), and the whole diffusion 10 pc region (Aharonian & Neronov 2005b). 
Interestingly, the angular scale of the TeV source was determined by HESS to be less than 
a few arc-minutes, indicating that this y-ray source is located in the central < 10 pc region 
(Aharonian et ah 2004). It suggests that the black hole Sgr A*, with a mass of 2.6 x IO^Mq 
(S chodel et al. 2002), should be involved (Aharonian & Neronov 2005a). 

Recently, a hadronic origin of TeV y-rays that is linked to the massive black hole has been 
addressed in detail by Aharonian & Neronov (2005b). They argued that the TeV y-rays are 
produced indirectly through the processes of 7r°-decay when relativistic protons are injected 
into the dense ambient gas environment. The flux of this radiation component depends on 
the density of the target, the diffusion speed of protons in the interstellar medium, and the 
injection rate of protons. As a result, a dense gas target and an extremely large proton flux 
is required. However, since the black hole Sgr A* currently only emits faint electromagnetic 
radiation, it is largely uncertain how the protons could be accelerated to relativistic speeds. 

It is well known that jets associated with accretion disks surrounding black holes are 
efficient in accelerating particles. For example, TeV y-rays observed from several BL Lac ob¬ 
jects (a subclass of active galactic nuclei) are argued to originate from relativistic jets (Pian 
et al. 1998). The jet model is also used to explain the production of TeV y-rays from micro¬ 
quasars (Bosch-Ramon et al. 2005). However, such a jet model cannot be applied directly to 
the black hole Sgr A*. An extraordinarily low bolometric luminosity of ~ lO^^er^fss”^ has 
been estimated for the black hole through multi-wavelength observations, which indicates 
that Sgr A* is in its quiescent dim state, and no powerful jet exists at present. 

We propose that the black hole Sgr A* could be re-activated and produce a powerful 
jet by capturing a red giant (RG) star, temporarily behaving like an active galactic nucleus 
(AGN). This may have happened thousands of years ago, naturally providing a mechanism 
to generate the proton flux required by the hadronic model of Aharonian & Neronov. The 
structure of our paper is as follows. The capture process is described in Sect.2. We estimate 
the model parameters and the energy of the proton flux available for the production of TeV 
y-rays in Sect.3. A brief discussion and conclusion are presented in Sect.4 
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2. Red Giant Capture and the Lighting up of Sgr A* 

2.1. the capture of a red giant 

Many if not all galaxies are expected to posses a massive black hole (Haehnelt & Rees 
1993). It has been thought that the capture and tidal disruption of stars by such a massive 
black hole could result in flarelike activities in the central regions of AGNs, and even in 
normal galaxies and globular clusters. The flare results from the rapid release of gravitational 
energy as the matter from the disrupted star plummets toward the black hole (Rees 1990; 
Shlossman et al. 1990). This process may even be the mechanism of y-ray bursts (Cheng & 
Lu 2001). It is interesting to note that the presence of an extended high-temperature plasma 
at the GC region has been conhrmed by Advanced Satellite for Cosmology and Astrophysics 
(ASCA) observations (Koyama et al. 1996), suggesting that the GC exhibited intermittent 
activities at least 300 yr ago. We propose that the black hole Sgr A* could have been 
naturally lit up and produced the currently observed TeV y-rays when it captured a RG 
star. 


The detailed capture and disruption process of a main-sequence (MS) star has been 
studied by several authors (Rees 1988; Cannizzo et al. 1990). Syer & Ulmer (1999) also 
discussed the capture rate of red giant stars by a massive galactic black hole in virialized 
star clusters. Compared with the capture of a MS star, the disruption of a RG star is different 
in that (1) RG disruption events generally last much longer, but are relatively fainter; ( 2 ) 
the tidal disruption radius of a RG star is typically equal to its pericenter radius; and (3) 
the size of the accretion disk formed in the capture process of a RG is larger. We believe 
that the capture of a RG star (rather than a MS star) is more appropriate to explain the 
observed y-ray emission from the GC (this problem is further addressed in section 4). 

We follow the treatments of Syer & Ulmer (1999) to study the capture of a RG star and 
scale our parameters to the physical conditions in the Galactic center. When a star of mass 
M* and radius R* passes by a black hole with a mass of Mbh, the star could be captured at 
an average tidal radius of 

Rt ~ ^ 8.35 X cm , (1) 

where m* = Mq = Mfe/i/ 10 ®MQ, and r* = R*/ 12 R 0 . For MS stars, R* = IR©, 

but for RG stars, the radius ranges from 3 to 2 OOR 0 (Syer & Ulmer 1999). The maximum 
radius of a RG star (200 R©) is attained only during a short period (10%) of its evolutionary 
lifetime, so the typical radius of RG stars is actually ~ I 2 R 0 , and we take r* = 1 in our 
model. The capture rate of MS stars by a massive black hole has been discussed by Phinney 
(1989), Rees (1990), and Cannizzo et al.(1990). However, the capture rate of RG stars 
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is more complicated, involving the mass function of stars, stellar evolution, the black hole 
mass, and the maximum radius of a RG (Syer & Ulmer 1999). By assuming a Salpeter mass 
function for the stars, Syer & Ulmer (1999) estimated the capture rate in our Galaxy as 
Nms ~ 4.78 X 10“^ yr~^ and Nug ~ 8.51 x 10“® yr~^ (Table 1 in Syer & Ulmer 1999) for 
MSs and RGs, respectively. 


Once a RG star is tidally disrupted, a fraction of its debris plummets to the black hole. 
The time for the debris to return to the pericenter (t^m) and the time for it to enter a 
circular orbit around the black hole {tdr ) are (Ulmer 1999) 


27iRl 


(GM,,)V2(2R,)3/2 


4.57m^ yr , 


( 2 ) 


tcir = norfetmm ~ 9.15(^^)(-^^^) yr , (3) 

2 A.blyr 

respectively, where Uorb is the number of orbits necessary for circularization, probably ranging 
between 2 and 10. We assume that the pericenter radius of the captured star equals its 
capture radius (Syer & Ulmer 1999) in the calculations. 

After the debris circularization, either a thick torus or a thin disk would be formed 
around the black hole, depending on two timescales. One is the time to radiate a signihcant 
portion of the energy of the bound debris at the Eddington limit {trad ), the other is the time 
to accrete a large fraction of the disk (face)- If Gcc ^ Uad, a thin disk will form; otherwise, a 
thick torus will appear (Ulmer 1999). The former case is necessary to our model, since the 
jet associated with the accretion disk is efficient to accelerate particles. Below we show that 
a thin disk will truly be formed in the process that we consider. 


The radiation timescale is (Syer & Ulmer 1999; Ulmer 1999) 

trad = ^ ~ 21^o.5?2r*Mg'^?7o.i yr , (4) 

L/^Ctet 

where ^ = 0.5.^o.5 is the fraction of the stellar mass that goes into the bounded debris (Rees 
1988), LEdd is the Eddington luminosity, and rj = O.lr^o.i is the efficiency of transferring the 
rest mass to radiation energy in the accretion process, which is dehned as ?7 = Ldisk/, 
where Ldisk is the total radiated disk luminosity, and M is the mass accretion rate of the 
disk. For a rapidly rotating black hole, rj ~ 0.42. 


For a standard accretion disk (Shakura & Sunyaev 1973), the accretion time can be 
estimated as. 


tacc = (^) 


,h. 


2.08 X 10 (-)_ 2 tt_ 3 m, / rJ yr , 




( 5 ) 
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where r and h are the radius and height of the disk, respectively, a is the viscous constant, 
and hlif = {GMhh /is the Keplerian velocity. The typical values for h/r are in the range 
10“^ to 10“^ (Shakura & Sunyaev 1973). We dehne {h/r )_2 = (h/r)/10“^ and q ;_3 = a/lO”^. 
From Eqs.(4) and (5), we hnd tacc ^ Uad, thus a thin disk forms in our model. 

The lifetime of our Galaxy is to ~ We hnd that when the black hole captures 

a MS or a RG star, the relation of maxitradRcirRacc) < min(l/iV, to) is satished, where 
N = Nms or Nrg- This implies that the activity of the black hole triggered by the capture 
events is intermittent. For the RG star capture, one long hare can be observed that lasts for 
tacc ~ 2.08 X 10^ yr. 


2.2. Accretion 


The time that the black hole becomes bright and active after the RG star capture 
depends on the evolution of its surrounding accretion disk. The overall evolution of the 
accretion rate is illustrated in Fig.l. Initially, the accretion rate of the disk exhibits an 
abrupt increase, since a local mass clumping in the inner-most orbit of the disk is similar 
to a delta-function (Rees 1998; see Fig.l). This stage will last until the mass accumulation 
rate of the disk reaches its peak Mpeak- We refer to this stage as phase 1. According to the 
simulations of Evans & Kochanek (1989), the value of Mp^ak is 

Mpeak ~ l.^Qrio.irnlr~^/‘^MY‘^^‘^MEdd , ( 6 ) 


where MEdd = 2.46 x 10 | Mq yr ^ is the Eddington accretion rate. 

When t > tpeak, the accretion rate evolves as (Rees 1988; Phinney 1989) 


1 M^ 
M ~ —- 


w5/3 


Q + . j- . ' 

^ ^mtn ^min 


2.97yo.imlr;^^^MY"^^MEdd{-^)-^^^ 

^m.in 


(7) 


From this point on, the black hole is activated and enters a luminous phase, acting like 
an AGN. This very active phase will continue until a critical time {tcrit) when the actual 
accretion rates falls below a critical value Merit (Beckert & Duschl 2004; 2002). We refer 
to this stage as phase 2. During this phase, the radiation efficiency of the accretion disk is 
very high. After phase 2, the accretion rate itself still continues to fall slowly as oc but 
the accretion flow becomes advection dominated. When advection takes over, the radiation 
efficiency drops drastically by several orders of magnitude (at least ~ 10^; Beckert & Duschl 
(2002)), resulting in a switch off of the AGN-like phase. We describe this stage as phase 
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3, during which the black hole becomes dim and no longer looks like a quasar. Beckert & 
Duschl (2002) pointed out that Merit should be less than 3 x because no consistent 

advection-dominated accretion flow (ADAF) models with M > 3 x are feasible. 

However, the exact value for Merit is still uncertain, since it is related to the radiation 
mechanism that is taking effect during the phase transition (Blandford & Begelman 1999; 
Quataert et ah 1999). In our study, we use the observed data to determine the value of Merit- 
Assuming that the transition from phase 2 to phase 3 is triggered by thermal instability of the 
disk (Shakura & Sunyaev 1973; Lightman 1974; Lu et al. 2000), the corresponding transition 
timescale is 


Attr ~ (aflx) ^ ~ 2.08a_3m^ yr . 


( 8 ) 


We see that the transition process is finished almost instantly. Note that a faint X-ray source 
(Baganoff et al. 2003) has been identified at the position of the GC, whose infrared radiation 
is also weak (Genzel et al. 2003). It indicates that the present accretion disk surrounding 
the black hole Sgr A* is probably ADAF (Narayan et al. 1995; Esin et al. 1998; Falcke & 
Markoff 2000). These observations also support our idea that the black hole is currently in 
phase 3. 

The time corresponding to Mpeak can be estimated from Eq.(7) as 

tpeak ^ ^■59tmin ~ 7.25m7 yr . (9) 


Although the jet formation itself is very difficult to model, observations are providing in¬ 
creasing evidence for the existence of powerful collimated outflows in black hole systems 
(Blandford & Konigl 1979; Fender et al. 2005; Gallo et al. 2005; Heinz & Grimm 2005). We 
assume that when the black hole Sgr A* evolves from phase 1 to phase 2, a powerful jet 
would form, originating at the inner edge of the disk (Markoff et al. 2001; Gheng & Lu 2001; 
Lu et al. 2000, 2003). For convenience, we define a dimensionless parameter 


Qj 


( 10 ) 


where Qj is the total jet power, including the rest energy of the expelled matter. The value 
of Qj can range from 6 x 10“^ to 0.4 (Fender et al. 2005; Gallo et al. 2005; Heinz & Grimm 
2005). 

It has been suggested that most of the black holes in the universe are rapidly spinning 
(Elvis et al. 2002; Volonteri et al. 2005). Recently, by analyzing the light curves of the X-ray 
and infrared flares from the GG region, the angular momentum of the black hole Sgr A* 
has been accurately determined as a = 0.99391 q;qq 74 (Aschenbach et al. 2004). The possible 
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mechanism for spinning up this black hole has been discussed by Y. Liu et ah (2006). 
Therefore, for Sgr A*, the radiation efficiency may be as high as 77 ~ 0.42. Combining 
Eqs.(7) and (10), we have 


t 


Qj = 1.75 X -) ^^^ergss ^ 

^min 

where 770.42 = 77 /O .42 and qj -2 = qj/lD~‘^. The total energy of the jet is then 

r^crit 


Ej^tot 


Qjdt 


^peak 


(11) 


( 12 ) 


The value of tcrit is determined in the next section. 


3. Jet Energy and TeV Gamma Ray Emission 

In spite of extensive observations, direct evidence for high-energy activities of the GC 
region, such as those observed in AGNs, is still lacking. Fortunately, strong fluorescent X- 
ray emission has been found from the cold iron atoms in the molecular cloud at the Sgr B2 
region (Sunyaev et ah 1993; Koyama et al. 1996). The most favorable explanation for this 
emission is that the molecular cloud was once illuminated by intense X-rays from Sgr A* 
(Sunyaev & Ghurazov 1998; Gramphorn & Sunyaev 2002). According to this theory, the 
X-ray luminosity of Sgr A* should have been Lx ~ 2 x lO^^ergs about 300 yr ago, which 
was most likely due to the capture of a star by the black hole. However, observations also 
show that the current X-ray luminosity of Sgr A* is only 2 x lO^^ergs (Koyama et al. 
1996). Therefore, if this explanation is correct, then the luminosity of Sgr A* has decreased 
rapidly by a factor of ~ 10^ in less than 300 yr. We note that in our framework, such a rapid 
dimming is possible only at the time of tcrit, when the accretion disk formed by capturing a 
star transitions from phase 2 to phase 3. As stated in Section 2 . 2 , the transition could be 
hnished in a few years (see Eq.( 8 )), and the luminosity could decrease by a factor of ~ 10^ 
after the transition. Thus, we suggest that the required luminosity of Lx ~ 2 x lO^^ergs 
was just the power of the black hole system when it was approaching the end-point of phase 
2. Then the lifetime of the jet can be determined by setting r]M{tcrit)c^ = Lx- The result is 

tcrit ~ 1-68 X 

Gombining Eqs.(9)and (11)-(13) and taking 77 = 0.42, qj = 10“^, we obtain 

Ej^tot — 2.76 X 10 ^^ ergs . ( 14 ) 



If a larger value of qj is assumed, then the Ej^tot would be even higher. Note that Ej^tot is the 
total kinetic energy of the jet. Possibly only ~ 10% of it could be converted into high-energy 
particles. 

The energy of the injected protons required by the hadronic origin of TeV 7 -rays has 
been modeled by Aharonian & Nerono (2005b). It depends on the diffusion coefficient of 
protons in the ambient gas. Generally, the dependence of the diffusion coefficient on proton 
energy can be expressed as D{E) = 10"^^{E/IGeV)^ k, 071 "^ where S ~ 0.3-0 .6 (Berezinskii 

et al. 1990) and k ~ 10“"^ to 10“^ (Aharonian & Atoyan 1996) are dimensionless parameters. 
There are mainly three typical propagation scenarios, characterized by different 6 and k 
values (Aharonian & Neronov 2005b). The case of the effective conhnement of protons (ECP) 
corresponds to 5 = 0.5 and n = 10“"^, the Kolmogorov type turbulence (KTT) corresponds 
to 5 = 0.3 and k = 0.15, and the Bohm diffusion (BD) corresponds to 5 ~ 1.0 and k ~ 10“^. 
For photons with a typical energy of 10 TeV, the diffusion coefficient is 1.0 x 10^® cm? s~^ 
for ECP, 2.4 X 10^® cm^ for KTT, and 1.0 x 10^° cm^ for BD. 

Assuming that the injection of protons lasts for ~ 10® yr, the required injection energy 
by the hadronic origin of the TeV 7 -rays in these three cases can be estimated as 

( 2.21 X 10^® ergs , for ECP , 

Ejtot = WplQ^yr ^ \ 2.38 X 10®°erc/s , for KTT , (15) 

\ 3.15 X 10®^ ergfs , for BD , 

where Wp is the required injection rate of protons corresponding to the three diffusion 
mechanisms, which is (Aharonian & Neronov 2005b), 


( 7.0 X 10®® ergss~^ , 

for ECP , 


Wp ^ < 7.5 X 10®’^ ergss~^ , 

for KTT , 

(16) 

\ 1.0 X 10®® ergss~^ , 

forBD . 



Comparing Eq.(15) with Eq.(14), and taking into account the efficiency (~ 10%) of convert¬ 
ing Ej^tot into cosmic rays, we hnd that the required energy could be satisfactorily provided 
by our RG star capture model in at least two diffusion cases, i.e., the ECP and KTT diffu¬ 
sion. In fact, if we assume qj > 0.1, then the jet energy would also be enough for the BD 
mechanism. 

Note that the linear size of the TeV source is likely to be only ~ 10 pc (Aharonian 
et al. 2004). We now further check which propagation scenario can meet this requirement. 
Theoretically, the protons would diffusion to a radius oi R = [4D(10TeV)Cjj]^/^ (Aharonian 
& Atoyan 1996), where tdif is the diffusion time of the injected protons when they propagate 
in the target. In our case, tdif = Grit + 300 pr. For the three diffusion mechanisms, the radius 
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of the source is, 


( 4.89 pc, for ECP , 

R = ^AD{10 TeV)tdif ^ < 75.9pc, for KTT , (17) 

y 489 pc , for BD . 

Eq.(17) shows that for both BD and KTT propagation scenarios, the corresponding radius 
is obviously too large as compared with observation. On the contrary, in the case of the 
ECP propagation scenario, we obtain R ~ 4.89pc, then the linear size is 2R = 9.78pc. It is 
in good agreement with the observed size of the TeV sonrce. We thus believe that the ECP 
diffusion is taking effect in our framework. 

In short, onr analysis shows that the captnre of a RG star by the black hole Sgr A* 
conld be the energy sonrce of the observed TeV y-ray emission from onr GC. The energy 
deposited in this way is large enongh to prodnce a relativistic ontflow of protons. The 
ontflow conld diffuse into a volnme with the diameter of ~ 10 pc. The detailed process that 
transfers the energy of the injected protons into that of observed y-rays has been discussed 
by other authors (Aharonian et ah 2004; Aharonian & Neronov 2005b), and is not addressed 
repeatedly here. 


4. Discussion and Conclusions 

A strong TeV y-ray source has been detected at the Galactic center, whose size is 
probably less than 10 pc. It should be closely related to the black hole Sgr A* of our Galaxy. 
Aharonian & Neronov (2005b) suggested a hadronic origin for these y-rays; i.e., they are 
prodnced throngh 7r°-decay process when a strong flow of relativistic protons interacts with 
the dense ambient gas. However, the nature of this proton flow is largely unknown. In this 
paper, we show that the reqnired proton ontflow could have been reasonably prodnced when 
the black hole Sgr A* captnred a RG star and formed an accretion disk aronnd it in the 
past. The whole process can be divided into three phases. In phase 2, which lasts for ~ 10^ 
yr, the black hole becomes active and Inminous thanks to a relatively high accretion rate, 
temporarily behaving like an AGN. A relativistic ontflow of protons can be ejected in this 
phase, whose total energy is as high as Ej^tot ~ 2.76 x 10®^ ergs, large enongh to meet the 
reqnirement of the Aharonian & Neronov model. 

According to our calculations, the propagation of protons in the target gas is through the 
EGP scenario. We show that the injected protons have diffusiond into a volnme of ~ 9.78 pc, 
consistent with the observationally inferred size of the TeV source. Furthermore, when the 10 
TeV protons diffnse in the target gas, their escape time is comparable with the characteristic 
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time of proton-proton interactions, implying that the spectrum of the observed TeV 7 -rays 
should be similar to that of the injected protons (Aharonian & Neronov 2005b). Since the 7 - 
ray spectrum observed by HESS is, J{E) = (2.5 ±0.21) x photon {cm? sTeV)~^ 

(Aharonian et ah 2004), we can conjecture that the initial spectrum of the injection protons 
should be Q{E) oc exp(—E/E q). Aharonian & Neronov (2005b) pointed out that the 
cutoff energy is Eq = 10 ^^ eV. 

In our study, we mainly consider the case in which the captured star is a RG star. If 
the captured star is a MS one, things will be much different. In Table 1, we have listed the 
key quantities calculated for a MS star capture, comparing them directly with those of a RG 
star capture. Interestingly enough, we hnd that the energy injection by a MS star capture 
also meets the requirement of energetics. However, we note that tcru > Gcc for a MS star; 
this means that the jet formed by the capture of a MS star terminates much earlier. In this 
situation, the diffusion timescale of the injected protons is tdif ~ face ± 300 yr, which is much 
smaller than the diffusion time of a RG capture event. Gonsequently, within the jet lifetime, 
the protons originating from a MS star capture will diffusion to 

( 1.05 pc, for the EGP , 

R = AD{lt)TeV)tdif ~ < 16.4pc , for the KTT , 

\ 105 pc, for the BD , 

which is inconsistent with the size (~ 10 pc) of the TeV source inferred from current obser¬ 
vations. This is the main reason that we prefer a RG star capture rather than a MS star 
capture in our framework. 

However, to determine the size of the TeV source observationally is not an easy task 
because of its extended nature. We note that the HESS collaboration has not yet published 
the hnal analysis of its observations of the Galactic center on larger scales, so there is still a 
lack of information on TeV emission beyond the 10 pc scale at the area. It is thus possible 
that the actual size of the TeV source may be larger. In addition, the reduction of TeV 
emission may also be due to the decrease of the density of the gas at larger distances, but 
not the lack of an ultra-relativistic proton flow. Taking into account these factors, the MS 
star capture scenario still cannot be completely excluded. In fact, a comprehensive analysis 
of the TeV radiation and a thorough investigation of the environment within ~ 100 pc around 
Sgr A* are necessary for us to understand the history of the activities in the GG region. 
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Table 1: Key quantities for a MS star capture and a RG star capture . 

tmin jyr) tpeakjyr) tacc jyr) tcrit {yr) Ejtotjergs) 


Star r.. 
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Fig. 1.— The behavior of the accretion rate versus time (M — t) when a black hole captures 
a star. The dashed line represents the local mass clumping in the inner region of the disk. 
The solid line is the evolution of the accretion rate as a function of 









